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In this paper, we study stochastic volatility models in regimes 
where the maturity is small, but large compared to the mean-reversion 
time of the stochastic volatility factor. The problem falls in the class 
of averaging/homogenization problems for nonlinear HJB-type equa- 
tions where the "fast variable" lives in a noncompact space. We de- 
velop a general argument based on viscosity solutions which we apply 
to the two regimes studied in the paper. We derive a large deviation 
principle, and we deduce asymptotic prices for out-of-the-money call 
and put options, and their corresponding implied volatilities. The re- 
sults of this paper generalize the ones obtained in Feng, Forde and 
Fouque [SI AM J. Financial Math. 1 (2010) 126-141] by a moment 
generating function computation in the particular case of the Heston 
model. 

1. Introduction. On one hand, the theory of large deviations has been 
recently applied to local and stochastic volatility models [1, 2, 4, 5, 20] and 
has given very interesting results on the behavior of implied volatilities near 
maturity. (An implied volatility is the volatility parameter needed in the 
Black-Scholes formula in order to match a call option price; it is common 
practice to quote prices in volatility through this transformation.) In the 
context of stochastic volatility models, the rate function involved in the 
large deviation estimates is given in terms of a distance function, which in 
general cannot be calculated in closed form. For particular models, such as 
the SABR model [19, 21], approximations obtained by expansion techniques 
have been proposed; see also [18, 22, 28]. Semi closed form expressions for 
short time implied volatilities have been obtained in [15]. 
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On the other hand, multi-factor stochastic volatility models have been 
studied during the last ten years by many authors (see, e.g., [8, 16, 18, 
27, 29]). They are quite efficient in capturing the main features of implied 
volatilities known as smiles and skews, but they are usually not simple to 
calibrate. In the presence of separated time scales, an asymptotic theory 
has been proposed in [16, 17]. It has the advantage of capturing the main 
effects of stochastic volatility through a small number of group parameters 
arising in the asymptotic. The fast time scale expansion is related to the er- 
godic property of the corresponding fast mean-reverting stochastic volatility 
factor. 

It is natural to try to combine these two modeling aspects and limiting 
results, by considering short maturity options computed with fast mean- 
reverting stochastic volatility models, in such a way that maturity is of 
order £<1, and the mean-reversion time, 5, of volatility is even smaller of 
order 5 = e 2 (fast mean-reversion) or 5 = e 4 (ultra-fast mean-reversion) . 

In [12], the authors studied the particular case of the Heston model in the 
regime 5 = e 2 by an explicit computation of the moment generating function 
of the stock price and its asymptotic analysis. 

In this paper, we establish a large deviation principle for general stochas- 
tic volatility models in the two regimes of fast and ulta-fast mean-reversion, 
and we derive asymptotic smiles/skews. For such general dynamics, a mo- 
ment generating function approach is no longer available. Our problem falls 
in the class of homogenization/averaging problems for nonlinear HJB-type 
equations where the "fast variable" lives in a noncompact space. We de- 
velop a general argument based on viscosity solutions which we apply to the 
two regimes studied in the paper. Viscosity solution techniques have been 
used in averaging of nonlinear H JB equations over noncompact space in [3] . 
However, the techniques in [3] were proved for a certain class of nonlinear 
HJB equations which does not include our case. In this paper, we develop 
a method more general than [3]. In particular, it can be used to treat the 
problems in [5], but not vice versa. 

We start by considering the following stochastic differential equations 
modeling the evolution of the stock price (St) under a risk-neutral pricing 
probability measure, and with a stochastic volatility determined by a pro- 
cess (Y t ): 

(1.1a) dSt = rStdt + a(Y t )S t dW t {1) , 

(1.1b) dY t = hm-Y t )dt + ^=YfdW t i2 \ 

o v 5 

where m € R, r, v > 0, and are standard Brownian motions with 

{W^\W^)t = pt, with \p\ < 1 constant. The process (Yt) is a fast mean- 
reverting process with rate of mean reversion 1/8 (5 > 0). The parameters /? 
and cr(y) are chosen to satisfy the following. 
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Assumption 1.1. We assume that: 




case of 1/2 < ft < 1, we require m > and Yq> a.s.; 
(3) a(y) e C(R;R+) satisfies 



a(y)<C(l + \yn 



for some constants C > and a with < a < 1 — ft. 

These assumptions ensure existence and uniqueness of a strong solu- 
tion of (1.1). This can be seen as a combination of existence of martin- 
gale problem solution (e.g., Theorem 5.3.10 in Ethier and Kurtz [9]) and 
the Yamada-Watanabe theory for 1-D diffusions (e.g., Chapter 5, Karatzas 
and Shreve [23]). In particular, Assumption 1.1(2) ensures that, in the case 
ft € [\, 1), Yt > a.s. for all i > (see Appendix A). In the case ft = 0, Y is 
an Ornstein-Uhlenbeck (OU) process with a natural state space (—00,00). 
In order to present both model cases using one simple set of notation, we 
denote the state space for Y as E$ with Eq := R if ft = and Eq := (0, 00) 
when ft £ [\, 1). 

Note that the Heston model, for which (3 = 1/2 and u{y) = y/y, does 
not satisfy these assumptions, but it has been treated separately in [12] by 
explicit computation of the moment generating function. 

The infinitesimal generator of the Y process, when 5 = 1, can be identified 
with the following differential operator on the class of smooth test functions 
vanishing off compact sets: 



Following the general theory of 1-D diffusion (e.g., Karlin and Taylor [24], 
page 221), we introduce the so called scale and speed measure of the (Yt) 
process, 



(1.2) 



B:={m-y)d y + \u 2 \y\ W dl 




— ^ ^- rl •y \ Tn ( 1 1 \ • — 




Denoting dS(y) := s(y)dy and dM(y) := m(y) dy, we then have 




Under Assumption 1.1 there exists a unique probability measure 



(1.4) 





m(y) dy < 00 



such that f Bfdiv = for all / G C^(E ). See Appendix C. 
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By a change of variable Xt = log St , we have 



dX t = (r - \a 2 (Y t )) dt + a{Y t ) dW\ 



(i) 



In order to study small time behavior of the system, we rescale time i i— >• et 
for < e <C 1; denoting the rescaled processes by X £t g j t and Y e g t , we have, 
in distribution, 



We are interested in understanding the two-scale e, 5 — > limit behavior 
of option prices and its implication to implied volatility. In this paper, we 
restrict our attention to the following two regimes: 



In view of [12], to obtain a large deviation estimate of option prices, it is 
sufficient to obtain a large deviation principle (LDP) for {X Et st'- e > 0}- By 
Bryc's inverse Varadhan lemma [7] (Theorem 4.4.2), we know that the key 
step is proving convergence of the following functionals: 



(1.6) u £ 4t,x,y):=elogE[e e ~ lh(x ^\X £AO = x,Y £A o = y], heC b (R), 



to some quantity independent of y. The rate function in the LDP is then 
given in terms of a variational formula involving the limit of the function- 
als u £:5 . 

For each h 6 C&(]R), the function u St s satisfies a nonlinear partial differ- 
ential equation given in (3.4). In Section 3.2, we use heuristic arguments to 
obtain PDEs that characterize the limit of these u e g- Proving this conver- 
gence rigorously, however, is nontrivial. Intuitively we know that, as Y has 
a mean reversion rate 1/5 and <5<S e, the effect of the Y process should get 
averaged out. To be exact, the form of nonlinear operator (3.5) indicates 
that convergence of u £ ^ is an averaging problem (over the fast y variable) 
for Hamilton-Jacobi equations. Such problems, in the context of compact 
state space for the averaging variable, can be handled by extending standard 
linear equation techniques using viscosity solution language. The Y process 
in this article lies in Eq, which is R in the case of (3 = and (0,oo) in other 
cases. Eq is a noncompact space, and therein lies an additional difficulty. 

We adapt methods developed in Feng and Kurtz [13]. Indeed, an abstract 
method for large deviation for sequence of Markov processes, based on con- 
vergence of H JB equation, is developed fully in [13] . The two schemes treated 
in this article are of the nature of Examples 1.8 and 1.9, introduced in Chap- 
ter 1, and proved in detail in Chapter 11 of [13]. In this article, we not only 



(1.5a) 



(1.5b) 




(i) 



5 = e 4 and 5 = e 2 . 
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present a direct proof, but also introduce some argument to further sim- 
plify [13] in the setting of multi-scale. This is possible in a large part due to 
the locally compact state space and mean-reverting nature of the process Y. 

In particular, modulo technical subtleties in verification of conditions, 
the setup of Section 11.6 in [13] corresponds to the large deviation result 
in our case of 5 = e 2 . Since Eq is locally compact, and we only deal with 
PDEs instead of abstract operator equations, great simplification of [13] 
can be achieved through the use of a special class of test functions. See 
Conditions 4.1 and 4.2. The techniques we introduce (Lemmas 4.1 and 4.2) 
are not limited to averaging problems, but are also applicable to problems 
of homogenization, which we will not delve into in this article. The rigorous 
justification of convergence of u Si s is shown in Section 5. 

The main results of the paper are stated in Section 2. Theorem 2.1 is 
a rare event large deviation- type estimate corresponding to short time, out- 
of-the-money option pricing. Corollary 2.1 and Theorem 2.2 give asymptotics 
of option price and implied volatility, respectively, for such situations. The 
proofs are given in the sections that follow, starting with heuristic proofs 
in Section 3.2 and finishing with rigorous justifications in Sections 4 and 5. 
The technical results in Lemmas 4.1 and 4.2 may be of independent interest. 

2. Main results. Observe that in the SDE (1.5), while the scaled log 
stock price process runs on a time scale of order e, the scaled Y process 
runs on a time scale of order e/S. This is due to the extremely short mean- 
reversion time, 5 = e r (r = 2,4), of the Y ei g. process. Thus, as e approaches 
zero, long-time behavior of the unsealed Y process comes into play. This 
long-time behavior of the Y process manifests itself in the large deviation 
principle (LDP) of the scaled log stock price via the quantities a 2 and H$ 
defined below. Define 

(2.1) a 2 := J a 2 (y)ir(dy); 

the average of the volatility function o 2 {-) with respect to the invariant 
distribution of Y . Recall B, the generator of the Y process, defined in (1.2). 
Define the perturbed generator 

(2.2) BPg(y)=Bg(y) + pavyPpd y g(y), geC 2 (E ). 
Let Y p be the process corresponding to generator B p , and define 

(2.3) H (p) :=limsup sup T' 1 lo gj B[e (1/2)|p|2 -fiT^O?)*^ = y]. 

Y p has strong enough ergodic properties that the limit above does not de- 
pend upon y even if we omitted the sup yg£ ; ; and, in fact, the lim supj.^^ 
can be replaced with ]imT-*x> in the above definition. We will justify this 
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fact in the rigorous derivations. By Girsanov's transformation 

(2.4) F (p) = limsupT- 1 log J B[ e /o T ^(^)^ (2) W+((iV)/2)H 2 /o T - 2 ^)^] ) 

T->+oo 

where Y is the process with generator B. From this expression, we see 
that Hq is convex and superlinear in p. H${p) is the scaled limit of the 
log moment generating function of a function of occupation measures of the 
process Y p . As such, it has an equivalent representation in terms of the 
rate function for the LDP of occupation measures of Y p . This equivalent 
representation of Hq is given in (5.12) in Section 5.2. 

Having defined these crucial terms, we proceed to the statement of our 
results. 

Theorem 2.1 (Large deviation). Assume X £i£ r^ = xq and Y £j£ t^ = y 
where r = 2,4 and suppose that Assumption 1.1 holds. For x G R, let 

I 1 2 

Xq — X\ 

(2.5) h(x;x ,t):= 2 _ 2 ^ , 
where a is defined in (2.1) and 

(2.6) I 2 (x;x ,t):=tL (^- 

where Lq is the Legendre transform of Hq defined in (2.3). 

Then, for each regime r € {2, 4}, for every fixed t > and xq £ R, pq € E$, 
a large deviation principle (LDP) holds for {X £t£ r t :e > 0} with speed 1/e 
and good rate function I r (x;xQ,t). In particular, 

(2.7) lim elogP(X £ e r t > x) = — I(x; xo,t) when x > xq. 

e— >0 ' ' 

Similarly, when x < xq, we have 

(2.8) lim e log P(X E , £ r tt <x) = -I(x; x ,t). 

Remark 2.1. The rate functions I r (x;xo,t), in both regimes, are con- 
vex, continuous functions of x and I r (xo;xo,t) = 0. 

Remark 2.2. In the case 5 = e 4 , observe that the rate function I^, 
in (2.5), is the same as the rate function for the Black-Scholes model with 
constant volatility a. In other words, in the ultra fast regime, to the leading 
order, it is the same as averaging first and then taking the short maturity 
limit. 

Remark 2.3. In the case S = e 2 , no explicit formula for the rate function 
is obtained. However, an explicit formula of the rate function is obtained 
for the Heston model in [12] which corroborates the formula in (2.6). The 
Heston model per se does not fall in the category of stochastic volatility 
models covered in this paper, but direct computation of Ho, given by (2.3) 
and Lq, its Legendre transform, is possible for this model. 
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Let So > be the initial value of stock price, and let X £t£ r = xo = log So- 
The asymptotic behavior of the price of out-of-the-money European call 
option with strike price K and short maturity time T = et is given in the 
following corollary. We only consider out-of-the-money call options by taking 

(2.9) S <K or x <logA". 

The other case, So > K, is easily deduced by considering out-of-the-money 
European put options and using put-call parity. 

Corollary 2.1 (Option price). For fixed t>0, 

lim elogE[e- r£t (S £ , £ r it - K) + ) = -I r (logK;x ,t) 

for r = 2,4. 

Denote the Black-Scholes implied volatility for out-of the-money Euro- 
pean call option, with strike price K, by cr r . )£ (t,log K, xq), where r = 2,4 
correspond to the two regimes. By the same argument used in [12], we get 
an asymptotic formula for implied volatility: 

Theorem 2.2 (Implied volatilities). 

(logK-x ) 2 



lim a r£ (t, log K,xq) 



6+ r ' £y ' & ' u/ 2I r QogK;x Q ,t)t 

Remark 2.4. In the case 5 = e 4 , the implied volatility is cf, which is 
obtained by averaging the volatility term o~ 2 (y) with respect to the equilib- 
rium measure for Y. It is likely that more features of the Y process, beyond 
its equilibrium, will be manifested in higher order terms of implied volatil- 
ity. Studying the next order term of implied volatility is a topic for future 
research. 

Remark 2.5. The limit of at-the- money implied volatility, that is, 
lim e _s.o of £ (t, xo, xo) , is obtained as in [12], Lemma 2.6. However, the conti- 
nuity of the limiting implied volatility at log K = xo is not obvious in the 
r = 2 case. We discuss this at the end of Section 6.3. 

3. Preliminaries. The process (X £ ^,Y £ ^) is Markovian, and can be iden- 
tified through a martingale problem given by generator 

A,sf(x,y) = e((r- ^cr 2 {y) J d x f(x,y) + ^a 2 (y)d 2 xx f{x,y) 

(3-1) 

+ ^Bf(x,y) + —f=pa(y)vyP d 2 xy f(x,y), 
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where / G C C 2 (R x E ). Recall that B is given by (1.2). Let g G C b (R) and 
define 

(3-2) v e)S (t, x, y) := E[g(X e>s>t ) \X e>Sfi = x, Y eA0 = y] . 

In gener al, v e> s G C b ([0,T] xlxE ). If, moreover, v £)S G C 1,2 ([0, TjxKxl), 
then it solves the following Cauchy problem in classical sense: 

(3.3a) d t v = A £>s v in (0, T] x E x £ ; 

(3.3b) v{0,x,y)=g(x), (x,y)eRxE . 

3.1. Logarithmic transformation method. Recall the definition of u £ ^ 
in (1.6). That is, u £)& ■= elogv £ j when g(x) = e £ ~ lh( - x \he C b (R), in (3.2). 
By (3.3) and some calculus, at least informally, (3.4) below is satisfied. This 
is the logarithmic transform method by Fleming and Sheu. See Chapters VI 
and VII in [14]. In general, in the absence of knowledge on smoothness of v £ ^, 
we can only conclude that u £ ^ solves the Cauchy problem (3.4) in the sense 
of viscosity solution (Definition 4.1). In addition to Fleming and Soner [14], 
such arguments can also be found in Section 5 of Feng [11]. 

Lemma 3.1. For h G C b (M), u £ ^ defined as in (1.6), is a bounded contin- 
uous function satisfying the following nonlinear Cauchy problem in viscosity 
solution sense: 

(3.4a) d t u = H £>s u in (0, T] xKx E ; 

(3.4b) u(0,x,y) = h(x), (x,y) G R x E . 

In the above, 

H £>5 u(t, x, y) = ee- e ~ lu A £)5 e £ ~ lu {t, x, y) 



(3.5) 



+ ha(y)d x u\ 2 + £ ^e~ £ ~ lu Be e - h 
2 o 



<vy p (A= d 2 xy u + d x u d y u 



where 



^re- £ ~ lu Be £ ~ lu = -Bu + 5~ l \\vy p 8 y u\ 2 . 
Note that H £ $ only operates on the spatial variables x and y. 



3.2. Heuristic expansion. By Bryc's inverse Varadhan lemma (e.g., The- 
orem 4.4.2 of [7]), we know that convergence of u e s is a necessary condition 
to obtain the LDP for {X £ ^ t : e > 0}. In this section, we describe heuristical- 
ly PDEs characterizing u £ s in the limit and the nature of convergence itself. 
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Henceforth, for notational simplicity, we will drop the subscript 5 and 
write u £ and H £ for u £i $ and H £ s, respectively. We begin by the following 
heuristic expansion of u £ in integer powers of e: 

(3.6) u £ = uq + eu\ + e 2 u 2 + e 3 U3 + e^u\ + ■ • ■ 

in both regimes. The Uj, i = 0, 1, . . . , are functions of t, x, y. In this heuristic 
section, we make reasonable choices of U{ which a posteriori, following a rig- 
orous proof of the convergence of u £ in Section 5, are shown to be the right 
choice. 

3.2.1. The case of 5 = e 4 . Computation of H £ u £ [see (3.5)] reveals that, 
in this scale, the fast process Y oscillates so fast that averaging occurs up 
to terms of order e 2 . Namely, no = uo(t,x), u\ = u\{t,x) and u 2 = u 2 (t,x) 
will not depend on y. To see this, we equate coefficients of powers of e in 
d t u £ = H £ u £ . 

Terms of 0{\) satisfy 

o = ^V /3 (^ ) 2 , 

so we choose uq independent of y. With this choice of u$ the equation for 
the coefficients of the next order terms, which is of O(Jj), reduces to 

= B Ul + \v 2 y w (d y ui) 2 . 

This equation is satisfied by choosing u\ independent of y. With this choice 
of u\, the equation for coefficients of the next order terms, of O(-), becomes 

= Bu 2 . 

By choosing u 2 independent of y the last equation is satisfied. 

Thus, by these choices of uq,ui and u 2 independent of y, it follows that 

H £ u £ (x,y) = \\a(y) d x u Q \ 2 + Bu 3 

+ e(a 2 {y) d x u d x ux + \o 2 (y) d xx u 

+ (r - \a 2 (y)) d x u + upa{y)y 15 d x u d y u 3 + Bu A ) 

+ o{e). 

The s° order terms then satisfy 

d t u (t,x) = \\d x UQ{t,x)\ 2 a 2 (y) + Bu 3 (t,x,y), 

that is, 

Bu 3 (t,x,y) = d t u (t,x) - l\d x u (t,x)\ 2 a 2 (y). 

The above is a Poisson equation for u 3 with respect to the operator B in 
the y variable. We impose the condition that the right-hand side is centered 
with respect to the invariant distribution ir [given in (1.4)]. This ensures 
a solution to the Poisson equation, which is unique up to a constant in y. 
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See Appendix B for growth estimates of the solution. Therefore we get 

d t u (t,x) = l\ad x u (t,x)\ 2 ; 

where 

a 2 = j a 2 (y)Tr(dy). 

Thus the leading order term in the heuristic expansion satisfies 
(3.7a) d t u = H u Q (x), t > 0; 

(3.7b) u (0,x) = h(x), 

where 

Hqu q (x) := \\a d x u G {x)\ 2 . 

3.2.2. The case of 5 = e 2 . When 5 goes to zero at a slower rate e 2 , limits 
become very different and more features in the Y process (rather than just 
its equilibrium) is retained. We observe that while no is independent of y as 
in the faster scaling regime, u\ may now depend on y. Equating coefficients 
of 0(e~ 2 ) in dtu £ = H £ u £ we get 

o=^y^ ) 2 , 

and so we choose uq = uo(t,x) independent of y. Then H £ u £ reduces to 
H £ u £ (t,x,y) = ~\a(y) d x u \ 2 + pa(y)vy r) ' d x u d y ui + e~ Ul Be Ul 

+ e(a 2 {y)d x u d x u 1 + \a 2 (y) d xx u + (r - \a 2 {y)) d x u 
+ Bu 2 + vy 2 ^ dyUi d y u 2 + pa(y)vy 13 d xy u\ 

+ po{y)vyP d x ui d y ui + pa(y)uy 15 d x u d y u 2 ) 

+ o(s). 

The leading order terms should satisfy 

d t u (t,x) = hd x u (t,x)\ 2 a 2 (y) + pva(y)y 13 d x u (t,x) d y ui(t,x,y) 

(3.8) 

+ e- Ul Be Ul (t,x,y). 

We will rewrite the above equation as an eigenvalue problem. Recall B, the 
generator of the Y process defined in (1.2) and the perturbed generator B p 
defined in (2.2). Then 

(3.9) e~ Ul Be Ul + pa{y)uy p d x u G d y ui = e'^ B 9xU °^ e u \ 

Fix t and x, and rewrite (3.8) in terms of the perturbed generator (3.9). 
e -«i B «.«o(t,*) e «i (t) X5 y) + l\ 9x u (t, x)\ 2 a 2 (y) = d t u (t, x). 
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Multiplying the above equation by e Ul , we get the eigenvalue problem 

(3.10) (B d ^ + V)g(y)=Xg(y), 

where V(-) = ^\d x uo(t, x)\ 2 a 2 (■) is a multiplicative potential operator, g(-) = 
e ui(t,x,-) anc j \(t^ x ^ _ d t u (t,x). Choose u\ such that (A, g) is the solution to 
the principal (positive) eigenvalue problem (3.10). Note that the dependence 
of the eigenvalue, A, on t and x is only through 8 x uq. If (3.10) can be solved 
with a nice g, then we have 

(3.11) X(t,x) = H (d x u ), 

where Hq is defined as (2.3). The leading order terms then satisfy 

(3.12) d t u (t,x) = H (d x u (t,x)). 

Constructing a classical solution for (3.10) is a considerably hard problem, 
even in the 1-D situation. If (3.10) can be solved with a nice g, then (2.3) 
always holds with the Hq given by (3.11). The converse is not always true. 
Especially, (2.3) says nothing about the eigenfunction g. However, we only 
need the definition in (2.3) in rigorous treatment of the problem. We will 
show (in Section 5.2) that (3.12) is the limit equation where Hq is given 
by (2.3) irrespective of whether a solution to the eigenvalue problem (3.10) 
exists or does not. 

To summarize, 

(3.13a) d t u (t,x) = H (d x u (t,x)), t > 0; 

(3.13b) u (0,x) = h(x), 

where Hq is given by (2.3) or (2.4). 

4. Convergence of HJB equations. The results of this section can be 
independently read from the rest of the article. 

We reformulate and simplify some techniques, regarding multi-scale con- 
vergence of HJB equations, introduced in [13]. Compared with [13], the sim- 
plification makes ideas more transparent and readily applicable. These are 
made possible because we are dealing with Euclidean state spaces which are 
locally compact. All these results are generalizations of Barles-Perthame's 
half-relaxed limit argument first introduced in single scale, compact state 
space setting. 

Let E C W n , Eq C W 1 and E' := E x E C R d where d = m + n. A typical 
element in E is denoted as x, and a typical element in E' is denoted as 
z = (x, y) with x £ E and y € Eq. We denote a class of compact sets in E' 

Q := {K x K: compact K CC E, compact K CC Eq}. 
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We specify a family of differential operators next. Let A be an index set and 

Hi(x,p,P;a):E x R m x M mxm xAh-E, i = 0, 1; 

H £ {z,p,P):E' xR d xM dxd ^R 

be continuous. For each / G C 2 {R d ), let Vf(x) G R d and Z> 2 /(aO € M dxd , 
respectively, denote gradient and Hessian matrix evaluated at x. We consider 
a sequence of differential operators 

H £ f(z) :=H £ (z,Vf(z),D 2 f(z)) 

for / belongs to the following two domains: 

D e j r :={/: / G C 2 (E'), f has compact finite level sets}; 

:= —D e)+ := { — f : f G C 2 (E'), f has compact finite level sets}. 

We will separately consider these two domains depending on the situation of 
sub- or super-solution. We also define domains D + ,D^ similarly replacing E' 
byE. 

We will give conditions where u e (t,z) = u £ (t,x,y) solving 

(4.1) d t u e (t, z) = H £ (z, Vu £ (t, z),D 2 u £ (t, z)) 
converging to u(t,x) which is a sub-solution to 

(4.2) d t u(t, x) < inf H (x, Vu(t, x),D 2 u(t, x);a) 

and a super-solution to 

(4.3) dtu(t, x) > sup H\(x, Vu(t, x),D 2 u(t, x); a). 

The meaning of sub- super-solutions is defined as follows (as, e.g., in Fleming 
and Soner [14]). 

Definition 4.1 (Viscosity sub- super-solutions). We call a bounded 
measurable function u a viscosity sub-solution to (4.2) [resp., super-solution 
to (4.3)], if u is upper semicontinuous (resp., lower semicontinuous), and for 
each 

u (t,x) = <j>(t) + f (x), 4> £ CHR + ), f €D + , 
and each xq G E satisfying u — uq has a local maximum [resp., each 

Ul {t,x) = 0(t) + h(x), 4> G C 1 ^), fx G L>_ , 
and each xo G E satisfying u — ux has a local minimum] at xo , we have 
d t uo(to,x ) - inf H (x , Vu (*o, ^o), D 2 u (t , x ); a) < 0, 

aSA 

respectively, 

d t ux(to,x ) - supH 1 (x ,Vu 1 (to,x ),D 2 ux(to,xo);a) > 0. 

agA 

If a function is both a sub- as well as a super-solution, then it is a solution. 
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We will assume the following two conditions. 

Condition 4.1 (limsup convergence of operators). For each /o G D+ 
and each a G A, there exists /o jE G Ar,+ (may depend on a) such that: 

(1) for each c > 0, there exists K x K £ Q satisfying 

{(x,y):H £ f 0j£ (x,y) > -c} n {(x, y) : / , e (x, y) < c} C K x A'; 

(2) for each A x A G Q, 

(4.4) lim sup |/o, e (a;,y)-/o(x)|=0; 

e ^ U (x,y)eKxK 

(3) whenever (x £ ,y £ ) £ K x K £ Q satisfies x e — > x, 

(4.5) limsupiA/ , £ (x £ ,y £ ) < H (x, V/o(x), D 2 f (x); a). 

Condition 4.2 (liminf convergence of operators). For each f\ G -D_ and 
each a G A, there exists /i )£ G A> e ,- (may depend on a) such that: 

(1) for each c > 0, there exists K x A G Q satisfying 

{{x,y) : H £ f h£ (x,y) < c} H {(x,y) : f 1:£ (x,y) > -c} C K x K; 

(2) for each A x A G Q, 

lim sup |/i(x)-/i |E (x,y)| = 0; 

(3) whenever (x £ ,y £ ) G A" x A G Q, and x e — >• x, 

limmf H £ f h£ (x £ ,y £ ) > H x {x,V f x {x),D 2 h{x);a). 

Let u £ be the viscosity solutions to (4.1); we define 
u 3 (t,x) := sup\limsupu £ (t £ ,x £ ,y £ ) :3(t £ ,x £ ,y £ ) G [0,T] x A x A, 

(t £ ,i £ )^(i,i),ii:xfeQ}, 

ii 4 (t,x) :=infhiminfu e (t e ,x e ,y e ):3(t £ ,2; e ,y e ) G [0,T] x A' x A", 

(t £1 i £ )4(t,4irx^eg}, 

and u = it| the upper semicontinuous regularization of 113 and u = (1x4)* the 
lower semicontinuous regularization of U4. 

Lemma 4.1. Suppose that sup £>0 || Me 1 1 00 < 00. Then: 

(1) under Condition 4-1, u is a sub-solution to (4-2); 

(2) under Condition 4-2, u is a super- solution to (4-3). 
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Proof. Let u (t,x) = <p(t) + f (x) for a fixed 4> G C 1 (]R + ) and / G D + . 
Let (to, xo) be a local maximum of u — uq, to > 0. We can modify /o and <p if 
necessary so that (io,£o) is a strict global maximum, for instance, by taking 
fo(x) = fo(x) + k\x — xo\ 4 and <fr(t) = (j)(t) + k\t — to\ 2 for k > large enough. 
Note that such modification has the property that 

lim sup |V/o(x)-V/o(x )| + |Z) 2 /o(x)-D 2 /oNI = 0. 

e ^°+|a;-xo|<£ 

Lettt o = + /o- 

Let q G A be given. We now take uo :£ (t,z) = 4>{t) + fo )£ (z) where /o, e is 
the approximate of /o in Condition 4.1. Since u £ is bounded, and UQ tE has 
compact level sets, there exists (t e ,z £ ) G [0,T] x E' such that 

(4.6) (u e -u 0:£ )(t £ ,z E )>(u E -u , £ )(t,z) for (t,z)e[0,T]xE' 
and 

(4.7) d t 4>{t e )-H E fo, E {z e )<Q. 

The above implies inf e H £ fo j£ (z £ ) > — oo. We verify next that fo, £ (z £ ) < 
c < oo. Then by Condition 4.1(1), there exists K x K G Q such that z e = 
(x E ,y E )eKxK. 

Take a (i, x) such that Cto(^) ^) < °°- Take z = (x, y) for some y G Eq. Then 
Uo,e(i, z) = 4>(i) + fo, £ (z) 4>(i) + fo(x) = u (i,x) < oo. 
Combined with (4.6), 

U Qi£ (t £ ,Z £ ) < 2 SUp 1 1 U £ 1 1 oo + SUp (M) < OO, 

£>0 £>0 

and sup e>0 fo, £ (z E ) < oo follows. 

Since K x K is compact in E', there exists a subsequence of {(t £ ,z £ )} (to 
simplify, we still use the e to index it) and a {Iq,xq) G [0, T] x S such that 
t £ — > to and x £ — > xq. Such (to,xo) has to be the unique global maximizer 
(to,xo) for u — uq that appeared earlier. This is because, by using x £ — > xq 
and z £ = (x £ ,y £ ), the definition of u and (4.4), from (4.6) we have 

(4.8) (u-uo)(to,x ) > (u-u )(t,x) V(t,x). 
Now, from (4.7) and (4.5), we also have 

d t u (t ,x ) < H (x ,V f (x ), D 2 f (x ); a). 

Note that to,xo and uq are all chosen prior to, and independent of, a. We 
can take inf Q£ A on both sides to get 

d t u (t ,x ) - inf H (x ,'Vuo(to,xo),D 2 uo(t ,x );a) < 0. 

The proof that u is a super-solution of (4.3) under Condition 4.2 follows 
similarly. □ 
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Lemma 4.2. Suppose that the conditions in Lemma 4-1 hold and that 
there exists h G Cb(E) such that 

lim sup \h(x) -u E (0,x,y)\ =0 Vif x K G Q. 

(x,y)£KxK 

Further suppose that for any sub-solution v>o(t,x) of (4-2) with uo(0,x) = 
h(x) and super- solution u\ of (4-3) with ui(0,x) = h(x), we have 

u (t,x) <ui(t,x), (t, x) G [0, T] x E. 

That is, a comparison principle holds for sub-solutions of (4-2) and super- 
solutions of (4-3) with initial data h. 
Then u = u = u and 

lim sup sup \u(t,x) — u e (t,x,y)\ = \/K x K G Q. 

5. Rigorous justification of expansions. To rigorously prove the conver- 
gence of operators H £ given by (3.5) to operators Hq obtained by heuristic 
arguments in Section 3.2, we rely on and extend results developed in [13]. 
An exposition of the relevant results from [13] was laid out in Section 4. 
In this section we verify Conditions 4.1 and 4.2 and prove the comparison 
principle in Lemma 4.2. We will adhere to the notation used in Section 4. 

Conditions 4.1 and 4.2 require us to carefully choose a class of perturbed 
test functions with an index set A and a family of operators {Hq(-; a), a); 
a G A} to obtain viscosity sub- and super-solution estimates of uq, the limit 
of u £ . This technique was first introduced in [13] and illustrated through ex- 
amples in Chapter 11 of that book. Our presentation simplifies the technique 
in the context of application here. We will make the sub-solution estimate 
given by Ho(-,a) tight, by inf-ing over a, hence introducing yet another 
operator Hq. Similarly, we sup over a to tighten up the super-solution type 
estimate provided by Hi(-,a) which introduces operator H\. 

Let 

(5.1) ((y) := |y-m| c , 

where C > is any number satisfying 2a < (, < 2(1 — /3) with a and /3 given 
as in Assumption 1.1. Throughout the two regimes (5 = e 4 ,e 2 ), we take the 
index set 

A:={a=(^6):^C 2 c (E ),0<e<l}; 
and define two domains 

D + := {/ : f(x) = <p(x) + 7 log(l + |x| 2 ); if G C 2 C (R), 7 > 0} 

and 

D_ := {/ : /(*) = <p(x) - 7 log(l + \x\ 2 ); <p G C 6 2 (M), 7 > 0}. 
A collection of compact sets in M. x Eq is defined by 

Q:={KxK: compact K CC R,K CC E }. 
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5.1. Case 5 = e 4 . For each / = f(x) G D + , and each a = (£,#) E A, we 

let 

g{y)--=t{y) + oay) 

and define perturbed test function 

f £ (x, y) := f(x) + e 3 g{y) = f(x) + e 3 £(y) + e 3 9((y). 
Note that \\dJWoo + ||dL/IU < oo. Then 

H £ f £ (x,y) = e[(r - \a\y)) d x f + \a\y) d 2 xx f] + \a 2 {y)\d x f\ 2 
+ B£(y) + 9B((y) + \e 2 v 2 y 2 ^\d y i{y) + 9d y ((y)\ 2 

+ ep<T(y)vyPd x f(d y t;(y) + d y ((y)). 

The choice of the number £ in definition of the function £(?/) in (5.1) guar- 
antees that BC,(y) < —C((y). Moreover, with the earlier assumption that 
< a < 1 — P, the growth of ((y) as \y\ — > oo dominates the growth in y of 
all other terms in H e f £ . Therefore, there exist constants cq,c\ > with 

HJ £ (x,y) < \\a{y)d x f{x)\ 2 + B£(y) - 9c ((y) + «*. 

In addition, 

f E (x,y) = f(x) + e 3 g(y) > f(x) - e 3 U\\oo- 

Furthermore, for each c > 0, we can find K x K £ Q, such that 

(5.2) {(*, y) : H £ f £ (x, y) > -c} D {(x, y) : f £ (x, y)<c}cKxK 

verifying Condition 4.1(1). The rest of Condition 4.1 can be verified by 
taking 

H (x,p;U) = sup ( l -\a{y)p\ 2 + B£(y) - 9c ((y) 

We define 

H f{x) := inf H (x,d x f(x);a) 

= inf inf sup (ha(y)d x f(x)\ 2 + B^y)-9c ay) 

Similarly, for / € a = (£,,9) E A, we can choose 
f £ (x,y) = f(x) + e 3 ay)-e 3 9((y). 
Then Condition 4.2 holds for the choice of 

ffiCs.p; £, 0) = inf ( \\a{y)p\ 2 + + flcoCfe)) • 
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We define 

Hif(x) := sa-pH 1 (x,d x f(x);a) 

aeA 

= sup sup inf ( ha{y)d x f(x)\ 2 + B£(y) + 9c ((y) 
o<e<i^c^(E )y^ E o V 2 

Next, to verify Lemma 4.2, we estimate H$f from above and H\f from 
below using some simple quantity. 

Lemma 5.1. 

H f(x) < \\ad x f{x)\\ f£D+; 

H 1 f{x)>\\ad x f{x)\\ /GL>_. 

We note that Hq,H\ have different domains D + and D_, respectively, D + (l 
D- =0. 

Proof. The key to obtaining the estimates in the statement of the 
lemma is the Poisson equation, 

(5-3) B X (y) = \\p\ 2 {a 2 -a\y)), 

where B is the differential operator (generator of Y) defined in (1.2). We 
will need growth estimates for \. In the case of f3 = (i.e., Y is an O-U 
process), Section 5.2.2 of Fouque, Papanicolaou and Sircar [16] contains 
such estimates. Specifically, if a(y) is bounded, \x(y)\ < C(l + log(l + |y|)); 
if cr(y) has polynomial growth, \ nas polynomial growth estimates of the 
same order. The following growth estimates for the situation | < j3 < 1 are 
derived in Appendix B: 

(5.4) |x'(y)l < Ci?/ 2cr_1 as y— > oo,for some positive constant C\. 

Therefore \x(y)\ < C(l + log(l + \y\)) if cr(y) is bounded and \x(y)\ < C*(l + 
y 2cr ) when < a < 1 — /3. 

We will make use of % as a test function in the expressions for HqJ 
and H\f . However, x does not have compact support. We choose a cut-off 
function (p to approximate it using localization arguments. Let nonnegative 
ip(y) £ C°°(Eq) be such that ip(y) = 1 when \y\ < 1 and when \y\ > 2. 
We take a sequence of £ n (y) = ^(^)x(2/)i which are truncated versions of x- 
Then 

Bt, n (y) = <p(j^jBx(y) + (m - y)x(y)n' 1 ip' 

+ \v 2 y w x{y)n- 2 v 




{I + ,VMy)"-V'(y. 
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Suppose a > 0. Noting that |y?(y)|, \(p'{y)\ and \<p"(y)\ are uniformly bounded 
and are when |y| > 2, and using the growth estimates (5.4) for \ an d x', 
we get 

i^(y)i <Wi+^+ (A 2 V+^WVV/^ 



< cy 2<J for all n. 

In the above, we used the fact that ^ < 2 and /3 — 1 < 0. Similarly, if cr(y) 
is bounded, that is, a = 0, we get |-B£ n (y)| is uniformly bounded for all n. 
Therefore, for large y, £(y) dominates B^ n {y) uniformly in n in the following 
sense: there exists a sub-linear function tp : M. \— > R + such that 

sup \BUy)\<^(((y))- 

n=l,2,... 

With the above estimate, we have 

H f{x) < limsup inf sup ( -\a(y) d x f(x)\ 2 + B£ n (y) - 9c ({y) 

n->oo 0<0<l ye E \Z 

<\\°d x f{x)\ 2 . 
Similarly, one can prove the case for H\f . □ 

By standard viscosity solution theory (e.g., [6]), the comparison principle 
holds for sub-solutions and super-solutions of 

d t u = \\ad x u \ 2 , t>0; 

u (0,x) = h(x), 

and the solution is uniquely given by the Lax formula (see [10]), 
(5.5) uo(t, x) = sup < h(x'" 



z'eR I 2a t 

Putting together the above result and Lemmas 4.1 and 4.2, we get: 

Lemma 5.2. 

lim sup \u £ (t,x,y) — uo(t,x)\ =0 Vc> 0, 
e ^°+|t|+M+|s/|<c 

where uq is the solution of (3.7) and is given by (5.5). 

5.2. Case 5 = e 2 . For each / = f(x) £ D + and a = (£, 9) £ A, we choose 
our perturbed test function as 

fe(x,y) := f(x) + eg(y), 
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where g(y) = (1 — #)£(y) + 0((y); ((y) is defined as before in (5.1). Then 
H e f E (x,y) = e[(r - ±a 2 (y)) d x f + \a\y) d 2 xx f] + ±e\y)\d x f\ 2 

+ e -9(y) B d*f(x) e g( y s ) 
< e[{r _ g xf + l a 2 {y) d 2 J} + l a l {y) \ dxf f 

+ (1 - 9)e^B d ^e^y) + 9e^B d ^e c (y), 

where B 9x -^ x ^ is the perturbed generator defined in (2.2). Recall that 
H^z/Hoo + ll^a;/lloo < oo by the choice of domain D + . We can thus find 
a constant cq > such that 

H e f e (x,y) < l\a(y)d x f(x)\ 2 + (1 - 9)e^B d ^e^y) + 9e^B dx -f^(y) + scq. 
Note that 

e-<B d ^e<(y) = Bay) + pa(yWd x f(x)d y ((y) + |^V^I^C(2/)| 2 , 
where 

(5.6) SC(y) = -C-\y- H c + ^ 2 y 2/3 C(C - l)\y - H c ~ 2 - 

The term — C(y) i n BC(y) dominates growth in y from all other terms in 
H £ f £ as |y| — > oo. Since £(y) — > oo as |y| — > oo, H £ f £ (x, y) — > — oo as |y| — > oo. 
We also have f £ (x,y) = f(x) + eg(y) > f(x) — e||£||oo- Therefore, for each 
c > 0, we can find K x K G Q, such that 

(5.7) {(x, y) : ff £ / £ (x, y) > -c} D {(x, y) : / e (x, y) < c} C AT x A' 

verifying Condition 4.1(1). 

The super-solution case follows similarly, where we define the perturbed 
test function as f £ (x,y) = f(x) +e(l + 0)£(y) — e9((y), for each / £ D_ and 

Take 

F o (x,p;e,0) := sup ( \\a{y)p\ 2 + (1 - e)e^B^(y) + for^e^y) ) , 

Htfap&e)- inf ('I| CJ (y)p| 2 + (l + ^) e -« J BV(y)-^ J BV(y)) 
yes \2 / 

and 

#o/0c) := inf inf H (x,d x f;£,9), 
o<e<i£eCg°(E ) 

H 1 f(x):= sup sup Hx(x,d x f;(,6). 

o<0<uec c °°(£ o ) 

Conditions 4.1 and 4.2 are satisfied by these choices of Hq and Hi. Note 
that, although ^|o"(y)p| 2 is not bounded in y, its growth is at most \y\ 2a and 
is dominated by the growth of C(y) for |y| large enough. 
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To verify Lemma 4.2, we develop useful sharp estimates for Hq and Hi 
next. Denote 

T(t)g(y) := E[g(Y t )\Y(0) =y], g £ C b (E ), 

and let B be the weak infinitesimal generator for semigroup {T(t) : t > 0} in 
Cb(Eo) (see page 244 of [13] for a definition of a weak infinitesimal genera- 
tor). Let D ++ (B) denote the domain of B with functions strictly bounded 
from below by a positive constant. Similarly define notations for B p , the 
weak infinitesimal generator corresponding to the process Y p introduced in 
Section 3.2.2. For each g G D ++ (W) C C b (E Q ), since C > 2<r, there exists 
compact K C C Eq with 

sup ( \\a{y)p\ 2 + (1 - 0) — (?/) + 0e-<B*>eHy)) 
yeE \ 2 9 J 

= sup ( ha(y)p\ 2 + (1 - 0) — GO + 0e-^^(y)) . 

For each e > 0, by truncating and mollifying g, we can find a £ := £ £ £ 
C™(E Q ) such that 

#o(z,p;£, 0) < e + sup f \\<r{y)p\ 2 + (1 - 0)^(y) + fle^S^y)') . 

Denote p = d x f(x). Then 

H f(x) < inf inf sup (\\a{y)p\ 2 + (1 - 9)^-(y) 

(5.8) 

+ 0e- c J B p e c (y)j. 

Similarly, we have 

Hif(x)> sup sup inf (I| (J (y)p| 2 + (l + 0)^(y) 
0<6»<l 5 GD++(BP)y e£; \^ 5 

(5.9) 

-Oe^BPe^ (y) J. 

We define Ib(sp) : ^(^o) ^ K U {+00} by 

I B U;p):=- inf / 5^d/zA / e-^ y) B p e^dn(y). 
geD++(MP)J Eo g J Eo 

However, we can find a sequence {g n } C D ++ (B P ) [take, e.g., y n := e^" where 
Cn £ C 2 (Eq) are some smooth truncations of £], such that 

/ e -^ y) B p e^ y) dfi{y) >limsup / 2^djt. 

J_B n— s-oo J£; 5Vi 
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Therefore we have 



(5.10) I B (ji;p) = - inf / —2-dfi. 

Recall that Y p denotes the process corresponding to generator B p (or, 
equivalently, M p ). It can be directly verified that Y p has a unique stationary 
distribution ir p and that Y p is reversible with respect to it (see Appendix C 
of this article). Let 



£ p (f,g):=-J 



fWgdn p 



be the Dirichlet form for Y p . By the material in Section 7 of Stroock [30] 
(particularly Theorem 7.44; note that the diffusion generated by B p has 
transition density with respect to Lebesgue measure, e.g., Theorem 4.3.5 of 
Knight [25]), we get 

,,n, lB(p;p) =e>(MM) = "4ryV ^ " 



dir p V dn p 



o 



see Appendix C.3 for the last equality above. If \i in Ib(h;p) is not absolutely 
continuous with respect to ir p , then the right-hand quantity in (5.11) is 
viewed as +oo. Again through Theorem 7.44 of [30], we also get that Hq, 
defined in (2.3), can be expressed as 

H (p)= sup ( J^j- / a 2 dfi- I B (fi;p)) 

(5.12) sup fM! f a 2 (y)h\y)7r p (dy) 

h£L2(7TP),\\h\\ L2(7rP) =l\ z JR+ 



u 2 



- j o y^\dh(y)\ 2 7r p (dy) 
As in Lemma 11.35 of [13], 

inf inf sup (h a (y)p\ 2 + (l-e)^-(y) + 9e^B p eHy)) =#o(p)- 

Using (5.8), this immediately gives 

H f(x)<H (df(x)), feD+. 
We will prove a similar inequality estimate for Hi , hence giving the following: 

Lemma 5.3. 

H 1 f(x)>H (df(x)), /€£>_, 
H f(x)<H (df(x)), f€D+. 
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It remains to prove the estimate for H\. By the proof of Lemma B.10 
of [13], 



(5.13) 



sup sup inf ( ha(y)p\ 2 + (l + e)^-(y)-ee-<BP e Hy)) 
0<6<i g eD++(np)y£ R +\ 2 9 J 

> inf lim inf t' 1 log E v [e^IPl 2 ti ds ] . 



We show that: 
Lemma 5.4. 

(5.14) limmit-hogE[e {1 ^ 2 ^ a2 ^ ds \YP = y]>H (p). 

Proof. The proof of (5.14) follows essentially the same argument used 
in Example B.14 in the Appendix of [13], which we will outline. Two in- 
gredients need to be emphasized. First, for each /i with Is(fJ,;p) < oo, by 
a mollification and truncation argument, we can find a sequence (i n (dy) = 
fe^n dwP d^iu) with h n + c n G C^°(Eq) for some constant c n , such that 
lim n ^oo I B (fJ> n ;p) = Ib(v,p)- Second, for every y £ E Q and every h £ C^°(E ), 
the following ergodic theorem holds: 



(5.15) 



where 



lim — E 

4— >oo t 



a 2 (Y s h )ds\Y h = y 



a\z) dn h (z), 



dY* = ((m - Y?) + ppa(Y*)v(Y s y + ^(W dh(Y?)) ds + V {Y?)» dWj, 

and where -fr h is the unique stationary distribution of Y h . We will prove (5.15) 
in Lemma 5.5. 

The process Y h is Y p under the Girsanov transformation of measures 



dP h 



dP 



■■expLpt) - h(Y ) - j\- h BPe h {Y s )ds\, 



where P and P h refer to the probability measures of the processes Y p 
and Y h , respectively. The invariant distribution of Y h is then 

e 2h dip 



dr 



We can write 
1 



lim inf — log E 

t— >+oo t 



f e 2h d<KV ' 
^\\p\ 2 f\ 2 {Y!)ds 



Y n p 



6 =V 
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1 n/i 

lim — log E 

t— >oo t 



ew<l\p\ 2 I <r 2 (Y s h )ds 



h(Y t h ) - h(Y h ) 



e- h B p e h {Y s h )ds 



y h =y 



1 „ph 

> lim -E p 

t— ¥00 t 



h(Y t h )-h(Y h ) 



e- h B p e h (Y s h )ds 



Yo=V 



(by Jensen's inequality) 
-\p\ 2 J j\z)d% h {z) 

/oo 
e~ h B p e h (z) dn h (z) (by ergodicity of Y ft ) 
-oo 



\\ P \ 2 j c °j\z)d* h {z)-e p \ 



dTr h / dfr' 1 



= -\p\ 2 / a 2 (z)^(z)-/(^;p). 

By arbitrariness of h, (5.14) follows. To complete the proof, we finally check 
that: 



Lemma 5.5. Equation (5.15) holds. 
PROOF. By Ito's formula, 

E[((Y t h )]=E[C(Y h )]+E 



B h C(Y s h ) 



ds, 



where B h ((y) = (m — y + ppa(y)vy 13 + v 2 y 2 ^ d y h(y))(' (y) + \v 2 y 2 ^"{y). As 
in (5.6), — C(y) is the dominating growth term in B h ((y). Therefore, defining 
a family of mean occupation measure, 



n h (t,y,A) :=E 



t~ 



1 {Y 3 h £A} ds 



Y h = y 
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we have that 
sup / ((z)ir h (t,y,dz) 

t>0 Jz 



■ supt" 

t>0 



l E 



C(Y s h )ds 











y 



<C(y;h(-))<oo. 



Hence {-k (t,y,-) :t > 0} is tight and along convergent subsequences and 
corresponding limiting point fr , we have 
ft 



(5.16) E 



Y 



h 



y 



(fd7T h , 



V G C b (E 



Such Tr h is necessarily a stationary distribution satisfying f B h tf)dTT h = 
for all i/j £ C 2 (Eq). Uniqueness of such probability measure can be proved 
by an argument similar to the one in Appendix C. We thus conclude that 
there is only one such Ti h and that convergence (5.16) occurs along the whole 
sequence, not just subsequences. Furthermore, the growth of a 2 is dominated 
by £, and so by uniform integrability argument, (5.15) holds. □ 

Now (5.9), (5.13) and (5.14) together give us the estimate for H\ in 
Lemma 5.3. 

From (2.4), we see that Hq{p) is convex in p e R. Let us denote its Leg- 
endre transform as Lq, then we have the following. 



Lemma 5.6. 
(5.17) 



The unique viscosity solution to (3.13) is 

— ( x — x' 



uo(t,x) := sup<J h(x ) 

x'el 



Moreover, u e converges uniformly over compact sets in [0,T] xMx£o to uq. 

PROOF. We know that no, defined by (5.17), solves (3.13) by the Lax 
formula. That uq is the unique solution follows from standard viscosity com- 
parison principle with convex Hamiltonians. The convergence result follows 
from multi-scale viscosity convergence results developed in Section 4, Lem- 
mas 4.1 and 4.2. □ 



6. Large deviation, asymptotic for option prices and implied volatilities. 

We finish the proof of Theorem 2.1, Corollary 2.1 and Theorem 2.2. 

6.1. A large deviation theorem. 

Proof of Theorem 2.1. From the previous section we have u E (t,x,y) — > 
uo(t,x) as e — > for each fixed (t,x,y) S [0,T] xlx Eq. All we need is ex- 
ponential tightness of {X Ej $ t t} to apply Bryc's lemma and to conclude our 
proof. This is obtained as follows. 

Let f(x) = log(l + x 2 ) and £(y) be defined as in (5.1). Take 

f c-r iA = I + e3 C(y)> for the case $ = e4 i 

{ ' \f(x)+eC(y), for the case 5 = e 2 . 
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Note that f(x) is an increasing function of \x\ and £(•) > 0; therefore, for 
any c > there exists a compact set K c Cl such that f £ (x,y) > c when 
x ^ K c . We next compute H £ f £ (x,y) [see (3.5)]. Observe that since ||5 a; /|| 00 + 
ll^xx/lloo < oo, by our choice of £(•), H £ f £ (x,y) — > — oo as \y\ — > oo. Therefore 
su ~Px&R,y&R H efe{x,y) = C < oo. For simplicity, we denote X eAt by X e>4 . 
The P and E 1 below denote probability and expectation conditioned on 
(X, Y) starting at (x,y). 

P(X £jt £ Kc ) e (c-fe(x,v)-tC)/a 

fe(Xe,t,Y e>t ) fe(x,y) 



< E 



cxp 



e e 

1 e -f e (X e , s ,Ye,s)/e Ase fe{X e , a ,Y e>a )/e dg 

Jo 

< 1. 

In the above inequalities, the term within expectation in the second line 
is a nonnegative local martingale (and hence a supermartingale) ; see [9], 
Lemma 4.3.2. We apply the optional sampling theorem to get the last in- 
equality above. Therefore 

e log P(X S)t i K c ) <tC + f £ (x, y)-c< const - c 

giving us exponential tightness of X e> t- 

Let iiQ ,r denote the limit of u e j when u Et s(0,x,y) = h(x) and 6 = £ r , 
r = 2, 4. Applying Bryc's lemma we get, {X £)£ r t } for r = 2, 4 satisfies a LDP 
with speed 1/e and rate function 

(6.1) I r (x;x ,t) := sup {h(x) - u^ ,r (t, x )}. 

h£C b (R) 

In Appendix D we check that l2(x; xq, t) = tLo( x °~ x ) where L is the Legen- 
dre transform of Hq defined in (2.3), and I4 = ^rpgr ■ O 

6.2. Option prices. 



Proof of Corollary 2.1. We follow the proof of Corollary 1.3 in [12] 
and show that lim e _ i> o+ £ ^°gE[(S £t t — K) + ] is bounded above and below by 
-I r (logK;x ,t). 

Recall that we are considering out-of-the-money call options and hence 
xq < log if [see (2.9)]. Since our rate functions I r (x;xo,t), for both r = 2,4, 
are nonnegative, convex functions with I r (xo;xo,t) = 0, they are conse- 
quently monotonically increasing functions of x when x > xq. Using this 
fact and the continuity of the rate functions, the proof of the lower bound 
follows verbatim from the proof in [12]. We refer the reader to [12] for details. 
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The upper bound follows from [12] once we justify the following limit: for 
any p > 1, 



(6.2) 



lim e\og E[S* S J = for both 5 = e 4 and 5 = e 2 . 



Recall the operator A e> s defined at the beginning of Section 3. By a slight 
abuse of notation, we can use A £ ^ to denote the operator acting on the 
unbounded function e px given below: 



A £ , S e px = e((r - \a 2 {y))pe px + \a 2 {y)p 2 e px ) . 



Let 



M t := expjpX M>t -pX £)5 , - J e~ pX ^ A £ ^e pX ^^ ds}. 

Then Mt is a nonnegative local martingale (supermartingale); this follows 
from the proof of [9], Lemma 4.3.2. By the optional sampling theorem, 



EM t < 1. 



Recall that X e g t t = logS £> s t t, then 



< (EM t ) 1/2 [E 



(6.3) 



exp|pA" £i5i0 + J 



e - P x eAs A P x sAs ds 



1/2 



(by Holder's inequality) 



< 1 • e pxa/2 E 



exp 



-pX s ,s,s Ae se PXe,S,s ds 







1/2 



We simplify and bound the right-hand side of the above inequality: 

ri 



E 



cxp 





E 

e erpt E 



e -pX EAsAe Se px eAs ds 

ex P |^^^r-^ 2 (y M)S )V+^a 2 (y M > 2 ') ds 
r ret/S 



by change of variable u = -s; recall that S = e or e' 

o 

et/8 



£rpt E expW-p) jT a 2 (Y u )dvj , 
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where Y u is the process with generator B given in (1.2). By convexity of 
exponential functions we get 



E 



(6.4) 



expj^ e- pX c^A £ , s e pX c 



< e £rpt E 



et/6 



exp{te(p 2 - p)a 2 (Y u )} du 



Since 5 = e 2 or e 4 , e/5 — > oo as e — > 0. Therefore, by the ergodicity of Y 
and exp{i(p 2 — p)<J 2 (y)} £ L l {dir) [this follows from an argument similar to 
proof of Lemma 5.5; note that a < 1 — f3 by Assumption 1.1(3)], the right- 
hand side of the above inequality (6.4) is uniformly bounded for all e > 0. 
Putting this together with (6.3), we get (6.2). □ 

6.3. Implied volatilities. 

Proof of Theorem 2.2. Recall that X £)t = logS £ ,t and xo = log So . 
Note that we have dropped the subscript 5 in the notation and the depen- 
dence on 5 = e 4 or e 2 should be understood by context. Our first step is to 
show that 

(6.5) lim a £ (t,logK, xo)Vei = 0. 

e^0+ 

Once we have shown this, the rest of the proof is identical to that of Corol- 
lary 1.4 in [12]. 

By the definition of implied volatility, 

log K + ret + a 2 et/2 s 



E[{S £ , t -K) + \ = e ret S Q <5> 



(6.6) 



Xq 



K<5> 



( 



Xq 



cr £ Vei 

log K + ret- a 2 et/2 
o- £ Vei 



where $ is the Gaussian cumulative distribution function. Let I > be the 
limit of a £ \fe~t along a converging subsequence. If lim £ _ s>0 + °f the left-hand 
side of (6.6) is 0, then I satisfies 



S $ 



Xq - log K l_ 



K J z -]agK _l_ 



I 



0. 



I ' 2, 

The only solution of the above equation is I = 0, and thus we get (6.5). 
We therefore need to prove 

(6.7) lim E[(S e , t - K)+] = 0. 

E-S>0+ 

By (1.5a) we have 

S £ , t -K = S -K + e I rS £ , t dt + y/i I S £ , t a{Ye,t) dW t {1) . 
Jo Jo 
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It can be verified that E[(S £! t - K) - (S - K)) 2 -)• 0, as e ->■ 0, for both 
cases 5 = e 4 and 5 = e 2 . Therefore 

lim E[(S e , t - K) + ] = E[(S - K)+] = 

as Sq < K (this is an out-of-the- money call option). 

The same formula is obtained when So > K by considering out-of-the- 
money put options. We finally turn our attention to at-the-money implied 
volatility. The asymptotic limit of at-the-money (ATM) volatility can be 
shown to be a 2 , that is, 

lim a 2 £ (t, log if, xq) = a 2 when xo = log A'; r = 2,4, 

£—7-0 ' 

by a similar argument as in [12], Lemma 2.6. The continuity, at-the-money, 
of the limiting implied volatility, that is, 

,. (logA-xp) 2 _ 2 
hm — — = a 

\\ogK-x \^Q 2I r (log K, x ,t)t 

is obvious in the r = 4 regime, but is more involved in the r = 2 regime. We 
conjecture that it is true, that is, 

z 2 

(6.8) lim — = = a 2 , 

and we briefly indicate an outline of the proof. Let 

A T (p) :=T- 1 logE[e/o T w CT (^)^i 2) +(( 1 -P 2 )/ 2 )N 2 /o T ^(^)^], 

so that Ho(p) = limy^oo A(p). The result (6.8) follows if Hq(p) is twice 
differentiate in a neighborhood of p = and Hq (0) = It can easily be 
checked that liniT'_ 5 . 00 A^(0) = The main difficulty is to get a uniform 
bound on A'^(p) for all T and in a neighborhood of p = 0. Obtaining such 
a uniform bound on Ay (p) involves tedious calculations but should follow 
from the multiplicative ergodic properties of the Y process (see [26]). □ 

In the following Appendix, we collect some material regarding 1-D diffu- 
sions Y and technical but elementary estimates. 

APPENDIX A: POSITIVITY OF THE Y PROCESS 

In this section we prove positivity of the Y process when i < /3 < 1 
in (1.1b). Assume m > and Yq> 0. Recall the scale function s(y) defined in 
the Introduction, and let S(y) = s(y) dy. By Lemma 6.1(h) in Karlin and 
Taylor [24], to prove that Yt remains positive a.s. for all t > 0, it is sufficient 
to show that 

lim S(e) = — oo. 
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For 0<e<l, 

-1 /•! 



/ / f f y 2(m — z) ~\ 

5 '( e ) = y s{y)dy = J expj-^ ^lup dz \ d V 

(where C is a positive constant and 2/3 — 1, 1 — @ > 0) 
l 

(positive integrand) dy 



2e 

f 2m 1 

- CgeXP I ^ (2 ^- 1)(2£ )2/3-l j^+ 00 

as e — >• + , provided m > 0. Therefore lim e _ s>0 + S(e) = —oo. 

APPENDIX B: GROWTH ESTIMATES FOR SOLUTIONS TO 

POISSON EQUATIONS 

Assume X satisfies the Poisson equation 

B X (y) = l\p\ 2 (a 2 -a 2 (y)), 

where a 2 , defined in (2.1), is the average of cr 2 (y) with respect to the invari- 
ant distribution ir(dy), given in (1.4), of the Y process. In this section we 
find growth estimates for x- 

The right-hand side of the above Poisson equation is centered with respect 
to the invariant distribution it{dy) = m ^ dy [given in (1.4)], and so 

roo 

(B.l) / m(z)(a 2 -a 2 (z))dz = 0, 

Jo 

where 

1 f fv 2(m - z) \ 

By (1.3), 

y 



X (y):= / dS(y) / \p\ 2 (a 2 -a 2 (z))dM(z 



o 
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is a solution up to a constant, and so 



x'(y) 



v 2 y 2 Pm{y) 

2 



IP I 



v 2 y 2 ^m{y) 



m(z)(a 2 — cr 2 (z)) dz 



m(z)(a 2 — cr 2 (z)) dz 



The last equality is by the centering condition (B.l). Given the bounds 
on <j(y) in Assumption 1.1(3), we can compute the following bounds where 
the constants, denoted by c, are positive and vary from line to line: 

12 



\x'(y)\ < 



c\p\ 



i/ 2 y 2 l 3 m(y) 
c\p\ 2 e 



z a m(z) dz 



2 p ay^ 



v 2 e -{y*-W)/{v*{l-P)) 
2m 



z 2a-2P e -az^ e -(z^)/(^(l-P)) 



where a = v i^~\\ > ^- Bounding e 



,1-2/3 



above by 1 we get 



\x\v)\ < 



clpPe * 1 - 2 " 



c\p\ 2 e' 



2a-2(3 e -(z 2 - 2 P)/(vi(l-P)) dz 



\2 p ay^ 



(by change of variable u = z 2 ~ 2/3 ) 



u (2«r-l)/(2-2^) exp 



U 



u 2 (l-P) 



du 



< 



elope " 1 " 2 ' 



y 2 *' 1 exp 



2-2/3 



^(l-/?) 



In the last inequality we used / a 00 [ tt ( 2ff - 1 )/( 2 - 2 « e - u /( ,/2 ( 1 -^)] du < i/ 2 (l 

(2 CT -l)/(2-2/3) p -a/(* 2 (l-/3)) 



/3)al 2a-l)/(2-2^, e -a/^(l- W) ( gince 2^1 <Q) _ Therefore 



lx'(s/)l < 



c|p| 2 e°* 1 " 9/ ' 



,2(7-1 



c|p| y 



2„,2<T-1 



as y — > oo, 



since e 



1 2/3 ~ 0(1) as y^oo. 



APPENDIX C: Y p PROCESS 

Fix pel. Denote // p (y) := (m — y) + ppa(y)uy 13 , and let Y p be the process 
with generator 

Bf S = ^(y) ^ 5 + \v 2 y^ d 2 yy g, g G C c 2 (£ ). 

In this section we calculate the unique stationary distribution and Dirichlet 
form of the process Y p , and we show that it is a reversible process. To this 
end, we first compute the scale function and speed measure. 
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The scale function and speed measure for the Y p process are given by 



s P (y) =exp|-^ 



2 7 2/3 



VZ 



and m p (y) 



Evaluating the integral in s p (y) we get (the C below denotes a positive 
constant that varies from line to line) 



s P (y) = < 



where 



C exp 
Cexp 



2m logy 

o T 



2-2/3 



2pp 



2m 



2-2/3 



+ 



y 



v 2 {2p-l)y 2 P~ 1 i/ 2 (l -/3) 



V 



if/3G0U ( -,1 ), 



J(y) 



o{z) 



dz. 



Due to bounds on a given in Assumption 1.1(3), there exist C\, C<i > such 
that 



where 

Therefore 
(C.l) 



C iy l -> 3 <J{y)<C 2y l -P + % 

0<1-/3<1-/3 + <x<1, if | <>S < 1, 
1 = 1-/3<1-/3 + cj<2, if = 0. 



— — — > when y — > or y — > oo, if - < /3 < 1, 
s P (y) 2 



k s p(y) 

Define for y e Eq, 



when |y| — ?• oo, 



SpCy) : = / Sp(2r)d2;. 



if p = 0. 



Observe that 5 p (y) — >• — oo as y approaches the left endpoint of Eq and 
<Sp(y) ~~ * +°° as y — » oo. 

C.l. Stationary distribution. Let 7r p be an invariant distribution of the 
process Y p . Suppose it has density function \P(y), that is, dit p {y) = ^?(y) dy, 
then \I/ is uniquely determined as the solution of 

l^ 2 y w m)--^MyMy)) = o 

satisfying *$>(y) > for all y and J^, $/(y)dy = 1. Solving the above differ- 
ential equation, we get ^(y) = m p (y)[CiS p (y) + C2]. Since ^ is nonnega- 
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tive, and S p (y) — > — oo as y approaches the left boundary of Eq, we take 
C\ = 0. The other constant C2 is uniquely determined by the condition 
J E ^{y)dy = 1. Therefore tt p is the unique invariant distribution of Y p and 
is given by 



(C.2) d^{y) = t>{y)dy 



Zi 



dy 



Ziv 2 y 2f3 s p (y) 



dy for y G E 



where Z\ = J Eq m p (y) dy. 

C.2. Reversibility. Let (p,tp € C 2 (E ), then 



Je Z\ J Eq 

=-( 

Z\ J En 



v 



5^^ + Mp(yy 



v 2 y 2 Ps p {y) 



dy 



<p e J 



d 



2,.2 / 3 



u-y 



dy 



dy. 



ZiJeq dy\s p (y)_ 
Integrating by parts twice and using the boundary conditions (C.l), we get 



i/jB p ipdir p 
Eq Zi 



Eq 



d v \ s p{y) 



dy 



pB P ^dTT P . 



E, 



C.3. Dirichlet form. By similar calculations as before, when proving re- 
versibility, we get, for f,g € L 2 (tt p ), 



£ p (f,g) 



fB p gd7r p 



Eq 



1 r m <(m)* 

ZiJeq dy\s p {y)J 



1 

,2 



Eq SpiV) 



dy 



-77 / y w f(y)g'{y)d^{y\ 



Eq 

where we integrated by parts once and used (C.l) in the second last line. 

APPENDIX D: RATE FUNCTION FORMULAS 
Recall the following characterization of the rate functions given in (6.1): 
I r (x;x ,t)= sup {h(x) - UQ ,r (t,x )}, 

h€C b (M) 



ASYMPTOTICS FOR STOCHASTIC VOLATILITY MODELS 33 

where r = 2, 4 correspond to the two regimes 5 = e 2 and 5 = e 4 , respectively. 
The Uq ,t are given in (5.17) and (5.5), respectively, as 

u (t, xo) = sup < n{x ) — tL\ 

hA,, \ r 7 / / \ / |^0 ~x'\ 2 



UQ A (t,x ) = sup \h(x') - (- 



2oH 



For notational convenience, we will drop the subscript r in I r , and, in 
Le case r = 4, we will denote the 
functions can then be rewritten as 



the case r = 4, we will denote the term ( ^ x °_^ ) by tL( x ° t x ). The rate 



I(x; xo, t) = sup inf < h(x) — h(x') + tL\ — ) > 

h€C h (R) x ' eR [ V t J) 



h&C b {R) 

for both regimes r = 2 and r = 4. 



Lemma D.l. 



I(x;x Q ,t) = tL^-j—\ 



Proof. Note that for both cases r = 2,4, Lq is convex, Lo(0) = and Lq 
is a nonnegative function. This is obvious for the case r = 4. We can deduce 
this in the r = 2 case since Hq(p) [defined in (2.3)] is convex and Hq(0) = 0. 

Re- write 



- / Xp-X 



I{x;x Q ,t) =tL \ — - 

+ sup inf Sh(x)-h(x')+tL (^^-) -tL (^^- 

where J = sup hgC . fc(R) J h and J h = mf x , m {h(x) - h(x') + tL (^^-) - t x 
Lo( x °~ :E )}. Taking x' = x in the inf we get J/j < and therefore 

(D.l) J<0. 

Note that xq and x are fixed. Define a function h* £ C&(M) as follows: 



hV)=t2o(^)AtLo(^ 



Then 
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and consequently 

(D.2) J>0. 
By (D.l) and (D.2), J = and we get 

I(x;x ,t) = tL ( 
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